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ABSTRACT 

Accreting black holes radiate in one of several spectral states, switching from one to 
another for reasons that are as yet not understood. Using the best studied example, 
Cyg X-l, we identify the geometry and physical conditions characterizing these states. 
In particular, we show that in the hard state most of the accretion energy is dissipated 
in a corona-like structure which fills the inner few tens of gravitational radii around 
the black hole and has Compton optical depth of order unity. In this state, an optically 
thick accretion disc extends out to greater distance, but penetrates only a short way 
into the coronal region. In the soft state, the optically thick disc moves inward and 
receives the majority of the dissipated energy, while the "corona" becomes optically 
thin and extends around much of the inner disc. The mass accretion rate in both states 
1 x 10- 8 Mt — - 1 



is 
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1 SPECTRAL STATES OF GBHC 

For many years we have known that Galactic Black Hole 
Candidates (GBHC) radiate X-rays in one of several spectral 
states, switching suddenly from one to another. Although in- 
termediate states are sometimes also seen, the spectral states 
of GBHCs are best typified by their extremes: a "hard" state 
(sometimes also called the "low" state because the 2-10 
keV flux is relatively weak in this state), and a "soft" state 
(sometimes also called the "high" state because the 2-10 
keV flux is relatively strong). The spectrum in both states 
may be roughly described as the sum of two components: a 
blackbody, and a power-law with an exponential cut-off. In 
the soft state the blackbody is relatively prominent and the 
power-law is steep; in the hard state more energy is carried 
in the power-law, whose slope is then shallower. 

More specifically, in the soft state (SS), the 0.1-5 keV 
band is dominated by a component best described as either 
a blackbody with temperature 0.3 - 1 keV, or perha ps a sum 
of blackbod ies all with temperatures in that range ( Mitsuda 
1984). The slope (in energy units) of the power-law 



et al. 



betw ee n 5 and 10 keV io generally 1.0 — 1.5 (Tanaka fe Low 



1995| ). OSSE and RXTE data, when fitted with a model 



consisting of a power-law with an exponential cut-off, show 
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slopes between 1.6 and 2.1 and a cut-off energy (although 
not very well determined) <; 200 keV (Grove et al. 1997, see 
also Phlips et al.1996, Cui et al. 1996 and Gierlihski et al. 
1997b for observations of Cyg X-l). Most of the energy is 
emitted in the soft thermal component. 

In the hard state (HS), the soft component is much 
less luminous, and can be re presented by a black body with 
a temperature 0.1-0 .2 keV ( Balucihska-Church et al. 1995 ; 
Ebisawa et al. 1996). Most of the energy is emitted in the 



hard tail, which can be represented by a power-law of slope 
0.3 - 0.7 with an exponential cut-off at about 100 keV 
(|Tanaka & Lewin 1995|; |Phlips et al. 1996|; |Grove et al. 1997|; 



Gierlihski et al. 1997a 



Zdziarski et al. 1997] ). The Gvnga 



spectrum of Cyg X-l also shows a clear sign of hardening at 
about 10 keV accompanied by a fluorescent Fe Kq line at ~ 
6.4 keV. These effects are most easily interpreted as the sig- 
nature of Compton reflection from cold matter in the vicin- 
ity of the hard X/7-ray source. In this interpretation, the 
strength of the Fe Kq line and the reflection bump requires 
the cold matter to subt end a solid angle £l/ 2n ~ 0.4 around 



the ha rd X-ray source ( Ebisawa et al. 1996 



1997a) 



Gierlihski et al 



It was originally thought that these systems were much 
more luminous in the soft state than in the hard. However, in 
a striking new result, Zhang et al. (1997b) have shown that 
the bolometric luminosity of Cyg X-l changes only slightly 
despite dramatic changes in spectral shape. In view of this, 
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we eschew the designations "high" and "low" , and strongly 
urge that in the future these states be called "soft" and 
"hard". 

What exactly is happening in these sources when they 
switch spectral state has long been a puzzle. However, in 
recent years advances have been made both on the observa- 
tional side (most importantly, hard X-ray/soft 7-ray obser- 
vations now give us a clearer picture of that portion of the 
spectrum) and on the theoretical side (we now understand 
much better how thermal Comptonization works when the 
seed photons are produced in large part by reprocessing part 
of the hard X-ray output). Making use of this progress, we 
can now use the observed character of the spectrum in these 
states to determine the geometry and energy dissipation dis- 
tribution in an accreting black hole system in all its different 
spectral states. The geometry and energy dissipation distri- 
bution inferred purely on the basis of radiation physics can 
then be used to guide efforts to obtain dynamical explana- 
tions for the changes in spectral state. Because the observa- 
tional data are best for Cyg X-l, we illustrate our method 
here by applying it to that source. As will be seen, even 
in this example, there are still uncertainties which prevent 
some of our conclusions from being as strong as might be 
possible in principle, but the bulk of this program is now re- 
alizable in the case of Cyg X-l, and it may soon be possible 
to extend it to other sources. 



2 DIRECTLY MEASURED PARAMETERS 
AND INFERENCES FROM ANALYTIC 
ARGUMENTS 

We attribute the two components with which the spectra 
are fit to two physically distinct, but related, regions: an op- 
tically thick, quasi-thermal region which is responsible for 
the blackbody component (the accretion disc?); and an op- 
tically thin very hot region which radiates the hard X-rays 
(a corona?). We call the intrinsic dissipation rates in the 
"disc" and "corona" Z/ S ntr and Lh, respectively. The scale 
over which the "disc" radiates most of its e nergy is R 



Table 1. Observational characteristics of Cyg X-l in its soft and 
hard states 



Following the consensus in the field ( ^ihajjuo^^light- 



man fc Eardlcy 1976|; |Haardt, Maraschi fc Ghisellini 1994} 



Stern et al. 1995; Zdziarski et al. 1997), we suppose that 



the hard X-rays are produced by thermal Comptonization. 
The seed photons for this process are partly created locally 
(by thermal bremsstrahlung or synchro-cyclotron radiation: 
Narayan & Yi 1995) and are partly created in the quasi- 
thermal region. The luminosity of the quasi-thermal region 
is in turn partly due to local energy dissipation, and partly 
due to reradiation of hard X-rays created in the "corona" . 

In this model, the shape of the Comptonized spectrum 
produced by the "corona" may be quite adequately de- 
scribed phenomenologically by two parameters: the power- 
law slope a, the exponential cut-off energy E c . Two more 
parameters (the effective temperature T s and L° bs ) define 
the soft part of the radiation. The ratio of the observed 
hard luminosity to the observed soft luminosity, L£ bs /I/° b8 , 
and the magnitude of the reflection bump C = Q/2tt (C — 1 
corresponds to the amplitude of reflection expected from a 
slab subtending a 2tt solid angle around an isotropic X-ray 
source atop the slab) complete the set of observables. These 
phenomenological parameters are determined by two dimen- 



Parameter 


Hard state 


Soft state 


^obs a 


1 X 10 37 erg/sl 1 ' 2 ' 3 ' 


4 X 10 37 erg/si 4 ' 5 ' 


robs b 
h 

a c 


4 X 10 37 erg/s 1 4 > 6 1 


1 X 10 37 erg/s ' 4 ' 5 ' 


0.6l 6 ) 


L6 [7,9,12] 


T s d 


0.13 keV' 1 * 2 ' 3 ] 


0.4 keVl 7 - 8 ' 9 ! 


E c e 


150 keV' 6 ' 10 - 11 ! 


> 200 keVl 7 ' 10 ! 


C f 


0.4(3,6] 


0.55' 12 ] 



a Observed soft luminosity (total luminosity below ~ 1 keV in the 

HS and below ~ 3 keV in the SS) 

b Observed hard luminosity 

c Energy spectral index 

d Temperature of the soft component 

°Cut-off energy of hard component 

'Covering factor of the cold matter 

References are:(l) Baluciriska & Hasinger 1991; (2) Balucifiska- 
Church et al. 1995; (3) Ebisawa et al. 1996; (4) Tanaka & Lewin 
1995; (5) Zhang et al. 1997b; (6) Gierlinski et al. 1997a; (7) Cui 
et al. 1996; (8) Belloni et al. 1996; (9) Dotani et al. 1996; (10) 
Phlips et al. 1996; (11) Zdziarski et al. 1997; (12) Gierlinski et al. 
1997b 

sional quantities, the total dissipation rate and R s , and four 
dimensionless physical parameters: the ratio L™ tT /Lh; the 
Compton optical depth of the "corona" tt; the fraction D of 
the light emitted by the thermal region which passes through 
the "corona"; and the ratio S of intrinsic seed photon pro- 
duction in the "corona" to the seed photon luminosity in- 
jected from outside. Another dimensionless parameter, the 
compactness lh = L^ot / '(rn e c 3 Rh), may be used to deter- 
mine the relative importance of pairs in the corona (Rh 
is the size of the corona). In this context it is also useful to 
distinguish the net lepton Compton optical depth r p from 
the total Compton optical depth (including pairs), tt. 

We will show how all these parameters, as well as several 
others of physical interest, may be inferred from observable 
quantities. The observables on which we base this analysis 
for Cyg X-l are shown in Table 1. Because Cyg X-l varies, 
the numbers seen at any particular time may be somewhat 
different from the "typical" values we cite. 

Certain of these figures deserve special comment. L° hB 
in the hard state is particularly uncertain. In works of 
Balucinska-Church et al. (1995) and Ebisawa et al. (1996), 
the spectrum in the ROSAT band is fit by a sum of a black- 
body and a power-law. They take L° bB to be only the part 
due to the blackbody; we believe that it is better described 
by the total of the two, and therefore find a soft luminosity 
two times greater than they do. In the soft state, C is very 
difficult to determine, since the amount of reflection depends 
on the assumed run of ionisation with radius (Chris Done, 
private communications). We will point out the degree to 
which our conclusions are sensitive to this uncertainty. 

Some of the physical parameters of the system may 
be derived (or at least constrained) almost directly from 
observables. For example, the electron temperature in the 
corona (measured in electron rest mass units) is very 
closely related to the cut-off energy of the hard compo- 
nent: O ~ f x E c /(m e c 2 ), where f x is a number slightly less 
than unity, whose exact value depends on tt. Similarly, 
L h » Lf s (1 + CA-Kda/dQ) /(l + C), where da/dQ, is the 
albedo per unit solid angle for Compton reflection, and we 
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D and C coefficients 




0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Figure 1. Coefficients D and C as function of h/Rh and R s /Rh- 
Solid curves represent contours of the constant D. Dashed curves 
give contours of constant C for Ty = 1.3 (corresponding to the 
HS). Dotted curves give C for tt = 0.3 corresponding to the SS. 
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Table 2. Inferred properties of Cyg X-l in its soft and hard states 



Parameter 


Hard state 


Soft state 


e a 


~0.2 


;> 0.3 


Lf^/L h b 


& 0.1 


~ 3 


Rs c 


500 km 


100 km 


D(l + S) d 


0.08 


0.3 


TT C 


~ 1.3 


^ 0.35 


R s /R h { 


~ 0.9 


~ 0.65 


h/R h z 


~ 0.2 


~ 0.5 


lH h 


~ 20 


~ 15 



a Electron temperature, B = kT e /mc 2 

b Ratio of intrinsic soft luminosity to the intrinsic hard luminosity 
(all the correction factors are here assumed to combine to unity) 
c Inner radius of the cold disc 

d Covering fraction of "corona" around thermal disc times the seed 
photon correction factor 

°Radial Thomson optical depth of the hot cloud 
Inner radius of the disc relative to the hot cloud 
g Elevation of the corona above to the disc relative to the size of 
the hot cloud 

h Hard compactness, lh = Lf l /Rh( a 'T/ m eC?) 



assumed that the hard component has an isotropic angular 
distribution and tt <C 1. 

The intrinsic disc luminosity is 



Lmtr rt 
e ' ' e 



CL h 
l + C 



1- / dttda/dn 



(1) 



where L\ ot = L° hs /p(6) [l - D (l - e~ TT )] is the total soft 
luminosity, and p(9) gives the angular distribution of the 
disc radiation (normalized so that J d£lp(8) = 4-k). Taking 
the disc emission to be approximately black body, its inner 
radius is 

1/2 

/.' - ! I ! ! . (2) 



\4ncjTt 



fr, 



where T 3 is the effective temperature at the inner edge. The 
additional correction factors are: / co i , the ratio between the 
local color temperature and the local effective temperature; 
/gr, which incorporates the general relativistic corrections 
linking the local color temperature to the observed one; and 
f m , which accounts for a proper integration over the disc's 
surface brightness distribution. For estimates of the correc- 
tion factors see, e.g., Shimura & Takahara (1995) and Zhang 
et al. (1997a). The combination of all these correction fac- 
tors is likely to be close to unity, but is uncertain at the 50% 
to factor of 2 level. 

We next employ the two following analytic scaling ap- 
proximations for thermal Comptonization spectra found by 
Pietrini & Krolik (1995): D(l + S) = 0.15« 4 L h /L* ot and 
tt = 0.16/(aO). The first expresses how the power-law 
hardens as the heating rate in the corona increases rela- 
tive to the seed photon luminosity; the second expresses the 
trade-off in cooling power between increasing optical depth 
and increasing temperature. 

Finally, we assume a description of the geometry in- 
volving the minimum number of free parameters consistent 
with describing a reasonable universe of possibilities. We 
take the disc to be an annulus of inner radius R 3 extending 
to infinite radius, with infinitesimal vertical thickness. We 
similarly assume that the corona is a pair of spheres uni- 



form inside, each of radius Rh, and centered a distance h 
from the disc midplane (when h = 0, the two collapse into a 
single sphere). This simplified description should be an ad- 
equate qualitative stand-in for any coronal geometry which 
is unitary (i.e. not broken into many pieces), axisymmetric, 
and reflection symmetric about the disc plane. Both C and 
D may then be written in terms of R„/Rh and h/Rh- We 
computed C taking into account electron scattering of the 
reflected radiation in the corona and assuming that the re- 
flection emissivity is tx (r 2 + h 2 )~ 3 ^ 2 (where Ti is the averaged 
height of the corona above the disc plane) and its angular 
distribution is oc cos#. In computations of D, we assumed 
that the disc surface brightness is oc r~ 3 and its intensity 
distribution is isotropic, i.e., p(6) = 2cos#. Contour plots of 
the coefficients C and D are presented in Figure ^. From C 
and D we can determine both R s /Rh and h/Rh. 

With these relationships and assumptions, we arrive at 
the derived quantities set out in Table 2. Although they are 
based on approximate expressions (and setting 4nda/dfl — 
0.3, along with p = /coi/gr//™. = 1)> we wm show in the 
following section that they are confirmed with only small 
quantitative corrections by detailed calculations. 

We can determine from Figure |i| that in the hard state 
h/Rh ~ 0.2. This corresponds to an almost spherical corona. 
The cold disc penetrates only a short way into the corona 
(see Gierlihski et al. 1997a and Dove et al. 1997 for a similar 
conclusion). If the disc were filled, relatively small C could 
be achieved for large r due to electron scattering of the re- 
flection component, but it would be difficult to get D small. 
The same argument also holds if the corona is larger than 
the disc (it would be easy to find geometries in which C is 
small, but then it becomes more difficult for D to also be 
small). In the soft state, h/Rh w 0.5, so that the corona 
is slightly elongated along the axis of the system (an out- 
flow?). Now, the corona covers a significant fraction of the 
inner disc (R s /Rh ~ 0.6). 

There is, of course, some uncertainty in determining C 
and D from observations. C is usually known with ~ 25% 
accuracy. This translates into similar uncertainty in R s /Rh, 
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and ~ 50% uncertainty in h/Rh- Estimations of D that are 
uncertain within a factor of 2 (mostly due to the strong 
dependence on the spectral index a) translate to ~ 30% and 
~ 15% uncertainty in R a /Rh for the SS and HS, respectively, 
while h/Rh is almost unaffected. 



3 DETAILED CALCULATIONS 

To verify these approximate arguments and refine our es- 
timates of the inferred system parameters, we have com- 
puted exact numerical equilibria for a variety of parame- 
ters. We solved energy balance and electron-positron pair 
balance equations coupled with the radiative transfer. This 
was done using the code described in Poutanen & Svens- 
son (1996). We assumed that the corona is spherical and 
centered on the black hole (i.e. h = 0), based on the fact 
that in both states the derived value of h/Rh is significantly 
smaller than 1. 

We divided these calculations into two sets, one de- 
signed to mimic the hard state, the other the soft. All but 
one of the system parameters needed to define an equilib- 
rium are given by the analytic estimates of the previous sec- 
tion. The one exception is t p , the net lepton optical depth. 
Although the total optical depth tt is reasonably well de- 
termined by observations, t v is not, because there is no easy 
way to tell from observations what portion of the Compton 
opacity is due to e pairs. In terms of the ratio z = n+ /n v , 
tt = t p (2z + 1). The one free parameter of the calculations 
is therefore t p . 

In the hard state models, we fixed lh = 20, T s = 
0.13 keV, Lf lr /L h = 0.1, and R s /R h = 0.9, and let t p = 0, 
1.0, and 2.0 (models HS1, HS2, HS3, respectively). Calcu- 
lations gave (tt, O, z)— (0.87, 0.39, oo) for model HS1, 
(1.12, 0.31, 0.06) for HS2, and (2.001, 0.16, 0.0002) for HS3. 
Selected predicted spectra are compared with observations 
in Figure ^. When t p <C 1, pairs dominate, and O ~ 0.4; 
T p 1 leads to 2 <C 1, and the greater optical depth de- 
presses the electron temperature. The best agreement with 
the observed spectrum (and the quality of this agreement 
is actually quite good) is found with tt — t v = 2, i.e., a 
state in which the pair contribution is negligible (assuming 
that pairs are thermal). Thus, the approximate parameters 
of Table 2 come very close both to a truly self-consistent 
equilibrium, and to producing an output spectrum like the 
one seen. 

In the soft state models we fixed lh = 15, T 3 = 0.4 keV, 
LT" /L h = 3, and R a /R h = 0.65, and let t p = 0, 0.1, and 0.3 
(models SSI, SS2, SS3, respectively). We find the following 
results: (r T , 6, z) = (0.24, 0.37, oo) for model SSI, (0.25, 
0.37, 0.75) for SS2, and (0.32, 0.31, 0.03) for SS3. 

In many respects, the soft state is quite similar to the 
hard state; it differs primarily in having roughly an order of 
magnitude smaller optical depth, and in the much greater 
importance of intrinsically generated soft photons. These 
changes do not alter very much the critical t v below which 
pairs dominate: it is still a few tenths. However, in contrast 
to the hard state, the best agreement with observations (and 
this fit is again quite satisfactory) is found with z ^ 1. 
Again we see that the analytic arguments are reasonably 
good guides to the results of more exact calculation. 



1 ^ 
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Figure 2. Theoretical spectra predicted by the calculations de- 
scribed in § 3. In both spectral states, the absolute normalization 
of the luminosity is arbitrary. The solid curve for the hard state 
corresponds to model HS3, and the dashed curve to model HS2 
(inclination of the disc is assumed to be 30°). The solid curve for 
the soft state corresponds to model SS2, and the dashed curve to 
model SS3. The discrepancies between predictions and observa- 
tions, particularly at the high energy end, for HS2 and SS3 show 
that there cannot be too many pairs in the hard state, or too few 
in the low state. The hard state data are from Gierlinski et al. 
(1997a) and the soft state data are publically available from the 
HEASARC. 



4 DISCUSSION 

These inferences allow several strong conclusions to be 
made. Most importantly, the fact that such a self-consistent 
picture may be drawn gives additional support to the basic 
picture of thermal Comptonization for the origin of the hard 
X-rays. This suggestion was part of the original proposal of 
Shapiro et al. (1976), who suggested that in the hard state, 
the inner portion of the accretion disc puffs up into a hot 
Comptonizing corona. Their picture of the geometry is also 
confirmed. In addition, we now see that in the soft state the 
disc extends inward almost as far as it can, while the corona 
surrounds it out to several gravitational radii. 

We are also now able to see that the disc receives only a 
minority of the dissipation in the hard state, but is the site 
of most of the heat release in the soft state. This conclusion 
is especially robust, for it depends primarily on the single 
ratio L£ bs /L° bs . If we nave overestimated L° hB in the hard 
state, the fraction of the dissipation taking place in the disc 
becomes even smaller. 

The radius of the disc's inner edge shrinks by roughly 
a factor of 5 between the hard state and the soft state. This 
ratio is somewhat sensitive to several uncertainties. First, 
R s cx [l/" 1 ^] 1 ^ 2 , and second, R s in the soft state can be 
affected by uncertainty in the scaling factors p(9), / co i, /gr, 
and f m . In terms of gravitational radii r g = GM/c 2 , R s ~ 
40r 9 M~ 3 1 in the hard state, and ~ 8r g M~^ in the soft state. 
Here Mio is the mass of the black hole scaled to lOM© , the 
best estimate of Herrero et al. (1995). 

Estimates of R s , combined with the intrinsic luminosity 
of the disc, also permit us to place strong bounds on the mass 
accretion rate through the disc, Md- Almost independent of 
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the spin of the black hole, mass passing through a thin disc 
extending outward from the ~ 40r 9 inferred for the hard 
state has a net binding energy ~ 0.02 of its rest mass. The 
intrinsic disc luminosity, at least, must be drawn from this 
store of energy; Md in the hard state must then be at least 
0.4 x 1O _S M0 yr _1 (this estimate is sensitive to uncertainty 

in L° s : it scales roug hly oc [L° bs ] ' ). The greatest possible 
Md would be required if all the energy for coronal heating is 
taken from matter accreting through the disc (that is, inside 
R s , the specific angular momentum of the matter changes, 
but not its specific energy). That upper bound is Md — 
4 x 10~ 8 Mq yr" 1 . In the soft state, the efficiency of the disc 
varies from ~ 0.05 (zero spin) to ~ 0.1 (maximum rotation). 
Similar reasoning then gives a lower bound of 0.7 - 1.4 x 
10 _8 Mq yr _1 and an upper bound only about 20% greater 
(because nearly all the luminosity in the soft state is radiated 
by the disc). If Md in the hard state is near the lower bound, 
the efficiency of the inner corona in the hard state must be 
~ 0.2 in order to explain the total luminosity. If so, the black 
hole would have to be rotating relatively rapidly. 

It is also of interest to compare Md to the accretion 
rate which would produce an Eddington luminosity if radi- 
ation were created at the efficiency of a maximal Kerr black 
hole — 6 x 10 -8 MiqMq yr -1 . In those units, the estimates 
of the previous paragraph give normalized accretion rates 
m = 0.07 (or possibly as high as m = 0.7) in the hard state, 
and ~ 0.1 - 0.2 in the soft state. 

Meanwhile, the corona shrinks in radial scale by only 
half as much as the disc when Cyg X-l moves from the 
hard state to the soft. In the latter state, it covers a sizable 
portion of the inner disc. Despite the sharp change in size 
and luminosity, the compactness of the corona changes by no 



more t han a factor of a few (incorporating the uncertainty 



ita 



in R s ) l ao a roault of the hard aoft tranoition, Howe 

optical depth drops by an order of magnitude. Although 
there are few e pairs in the equilibrium characterizing the 
hard state, they can be comparable in number to the net 
electrons in the soft state. This fact also means that the fall 
in t p from the hard state to the soft is even greater than the 
fall in tt- 

We close by noting that in the hard state there must be 
a hole in the central region of the thermal disc which is much 
bigger than the radius of the innermost stable orbit. The ori- 
gin of this hole is a question of prime importance. One possi- 
bility arises from the fact that the disc becomes dominated 
by radiation pressure inside ~ 35(m/0.1) 2,/3 gravitational 
radii. As shown by Shakura & Sunyaev (1976), in the strict 
Q-model, radiation pressure-dominated discs are thermally 
unstable. It is interesting that the inner radius of the disc 
in the hard state corresponds quite closely to this possible 
point of instability. We are then left with the question of 
how a thin disc is able to persist in to much smaller radius 
in the soft state, when m is still ~ 0.1. Another possibility is 
that the inner disc in the hard state is advection-dominated, 
while in the soft sta te the accretion rate is just a bit too high 
for that to happen (Narayan & Yi 1995). While the creation 



of an advection-dominated region may be a nonlinear end- 
state of the thermal instability, so that this model provides 
an attractive explanation for the hard state, it also leaves 
open how the thin disc survives thermal instability in the 
soft state. 
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